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ABSTRACT 

— 

Numerical  solutions  to  the  equations  which  govern  surface  water 
waves,  generated  by  axial ly-symmetric  initial  impulses  applied  to 
the  surface  of  a water  medium,  are  presented  in  this  paper.  These 
equations  were  first  derived  by  Kranzer  and  Keller  in  their  paper, 
"Water  Waves  Produced  by  Explosions." 

Several  cases  are  considered  within  the  paper.  The  initial 
impulse,  depth  cf  medium,  and  time  and  distance  from  the  impulse  are 
varied.  Finally,  comparison  is  made  between  the  theoretically  pre- 
dicted waves  and  the  waves  measured  at  two  different  test  conditions. 
In  both  cases,  fairly  close  agreement  was  found  between  the  theoreti- 
cal and  empirical  data. 
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Nomenclature 


A variable  portion  of  the  amplitude  of  the  wave  envelope 
Amax  maximum  value  of  A 

p 

g gravity  constant  (taken  as  .00980  Km/ sec  ) 

h uniform  finite  depth  of  the  medium 

p 

I0  value  of  the  initial  impulse  at  the  origin  (dyne  sec/cm  ) 

T(o/h)  zero  order  Hankel  transform  of  the  initial  impulse  distribution 
Jx(y)  the  xth  order  Bessel  function  of  y 
Q energy  content  of  the  initial  impulse 

R effective  radius  of  the  initial  impulse  (kilometers) 

r radius  from  the  origin  of  the  impulse  (kilometers) 

T period  of  the  wave  (seconds) 

t time  measured  from  the  initial  impulse  (seconds) 

n(r,t)  height  of  the  surface  wave  at  distance  r from  the  origin, 
and  at  time  t 

A wavelength  of  the  wave  (kilometers) 

O 

p density  of  the  medium  (taken  as  1.025  Gr/cm  ) 

o dimensionless  auxiliary  variable 

♦ dimensionless  auxiliary  variable 

4>'  derivative  of  $ with  respect  to  o 
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INTRODUCTION 

In  this  paper  numerical  solutions  are  given  for  the  equations  which 
govern  water  waves  generated  by  an  axial ly-symmetric  initial  impulse  which 
is  applied  on  the  surface  of  a water  medium.  These  equations,  developed 
by  Kranzer  and  Keller,1  are  valid  in  the  far  field  and  for  a medium  of 
finite  depth  h.  The  waves  generated  will  obey  the  linear  theory  of 
surface  waves. 

The  initial  impulse  imparted  to  the  uniform  medium  is  axial ly-symmetric, 
and  the  resulting  axial ly-symmetric  surface  water  waves  emanate  from  the 
source  radially.  The  height  of  the  water  wave  is  denoted  n(r,t),  where 
r is  the  distance  from  the  origin,  and  t is  the  time  from  the  initial 
impulse.  The  equations  developed  by  Kranzer  and  Keller  involve  an  asymp- 
totic expansion  and  use  the  Method  of  Stationary  Phase  (first  applied  to 
this  problem  by  Lord  Kelvin  in  1887),  valid  for  large  values  of  r and  t, 
and  can  therefore  only  be  applied  in  the  far  field,  i.e.,  for  r » R, 
where  R is  the  effective  radius  of  the  impulse. 

Kranzer  and  Keller  assumed  that,  "the  initial  impulse  at  any  point  in 
the  surface  is  determined  by  the  impulse  in  the  incident  shock  wave  at  the 
point."  It  was  also  assumed  that  the  initial  displacement  of  the  entire 
medium  is  zero.  In  this  analysis  only  the  case  of  a surface  impulse  is 
considered. 

A computer  program  was  developed  for  a CDC  7600  computer  located  at 
the  Lawrence  Berkeley  Laboratory.  With  either  a specified  value  of  r or 
t,  the  program  generates  values  of  n(r,t).  The  Input  parameters  Include 
oT^the  medium  (assumed  to  be  sea  water),  depth  of  the  medium,  radius 
of  the  charge,  and  Initial  impulse  distribution.  A graphical  display  of 
these  data  is  possible  using  the  6DS  (Graphical  Display  System)  available 


1 


2 

in  the  computer  library.  Several  postprocessors  are  available  for  the  GDS, 
and  the  graphs  in  this  report  are  of  the  microfiche  type. 

DISCUSSION  OF  THE  KRANZER  AND  KELLER  EQUATIONS 

The  medium  under  consideration  is  one  consisting  of  an  incompressible 
fluid  which  has  an  upper  free  surface  and  a lower  rigid  surface  at  a constant 
depth  h.  The  boundary  of  the  medium  is  taken  at  infinity.  This  implies 
that  any  waves  produced  will  be  a direct  effect  of  disturbances  from  within 
the  medium,  and  not  from  any  reflections  from  the  boundaries.  For  the 
source  of  a disturbance,  a cylindrically  symmetric  distribution  of  impulse 
imparted  at  the  origin  (i.e.,  r = 0)  at  t = 0 will  be  considered.  The 
result  of  this  impulse  will  be  the  generation  of  cylindrically  symmetric 
surface  waves  emanating  from  the  origin.  It  is  assumed  that  the  upper 
surface  has  no  initial  displacement  nor  initial  velocity.  At  any  time  t 
and  any  distance  r for  the  origin,  the  wave  height  is  described  by  the 
function  n(r,t). 

Beginning  with  the  Laplace  equation  which  satisfies  the  potential 
function  of  the  flow,  and  applying  the  Hankel  transform  to  this  equation, 
then  according  to  a linear  theory  of  waves,  n(r,t)  is  given  as: 

n(r,t)  * i -jx  lim  / S3^2  T(s)(tanh  sh)^2  sech  sh 

pg  ' y -*•  o’  o 

x cosh  (y+h)  sin  [(gs  tanh  sh)1^2  t]  JQ( rs)  ds  (1) 

where  s is  defined  as  o/h.  For  a derivation  of  this  equation,  see 
Stoker.2 

In  Eq.  (1)  the  function  T(s)  is  defined: 
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I(s)  = / I (r ) J (sr  r dr  (2) 

o 

The  integral  in  Eq.  (1)  is  evaluated  by  use  of  an  asymptotic  expansion. 
It  is  also  necessary  to  expand  the  Bessel  function  by  its  asymptotic 
expansion  valid  for  large  arguments;  i.e.. 


J0(*>  = ( if  )1/2  {Po^x^  cos  (x  - J ) - Q0(x)  sin  (x  - J )} 


(3) 


where 


,2  ,2  ,2  -2  ,2  72  ,2  ,2  ,2  72  02  ,,2 

p ~ -j  1 »3  + 1 «3  *5  •/  _ 1 *3  *5  *7  «9  •!!  + 

O Ol/O..  \ 2 * i / o..  Cl/O..  \6 


2! (8x)d  4! (8x) 


6!  (8x) 


l2  + 12»32,52  _ 12»32>52»72»92  + 
l!^8x^  3 ! ( 8x ) 3 5 ! ( 8x ) 5 
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Then  it  is  possible  to  use  Kelvin's  Stationary  phase  formula  to 
arrive  at  Eqs.  (4)  to  (8).  In  this  case,  the  argument  is  (rs)  and  the 
expansion  is  valid  when  r is  large,  i.e.,  only  in  the  far  field. 

Instead  of  using  two  variables,  it  is  assumed  that  the  ratio  r/t  is 
fixed,  and  then  only  one  large  parameter  will  exist,  either  r or  t. 

The  arguments  of  the  trigonometric  functions  for  the  expansion  of  the 
Bessel  function  contain  one  term  which  is  very  large,  r,  and  one  which 
varies  more  slowly,  s = o/h.  So  as  s varies  slightly,  the  trigonometric 
function  will  go  through  many  oscillations.  Assuming  that  the  coefficient 
in  front  of  the  trigonometric  functions  is  somewhat  bounded,  then  the 
positive  areas  under  the  cosine  curve  will  be  cancelled  by  the  negative 
areas  when  the  integral  is  taken  from  zero  to  infinity.  However,  there 
may  be  some  places  where  the  oscillation  is  much  slower,  and  the  cancella- 
tion will  not  be  complete.  These  points  can  be  shown  to  occur  when  a zero 
of  one  of  the  derivatives  of  s occurs  within  the  limits  of  integration. 
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s is  expanded  about  that  point,  and  the  coefficient  of  the  trigonometric 
function  is  evaluated  at  the  point.  Another  important  element  to  be  noted 
is  that  the  sine  is  rewritten  in  terms  of  the  cosine,  and  it  is  possible 
to  write 


and  a Itself  is  a dimensionless  auxiliary  variable.  The  sign  'v*  has  been 
used  to  indicate  "asymptotically  equal  to." 

In  Eqs.  (4)  to  (8),  IQ  Is  the  Initial  Impulse,  R is  the  effective 
radius  of  the  Impulse,  p is  the  density,  g is  the  gravity  constant. 
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r is  the  distance  from  the  origin,  t is  the  time,  A ts  the  variable 
portion  of  the  amplitude  of  the  wave  envelope,  T is  the  period,  X is 
the  wave  length,  h is  the  depth  of  the  medium,  and  T is  given  by 
Eq.  (2). 

It  should  be  noted  that  Eq.  (7)  is  Eq.  (2.5)  of  the  Kranzer  and  Keller 
paper  except  for  the  o in  Eq.  (7).  This  sigma  is  consistent  with  Kranzer 
and  Keller's  Eq.  (2.11)  and  also  agrees  with  Eq.  (B.l)  of  a paper  written 

4 

by  Kriebel.  There  exist  other  errors  in  the  original  Kranzer  and  Keller 
paper  which  were  suspected  and  later  verified  by  the  Kriebel  paper.  The 
graphs  which  are  given  as  Figs.  1,  2,  and  3 display  |A|,  <j>(o),  and 
[tanh  a <p(a)/~  <p'(o)]  ' . These  diagrams  are  taken  from  the  Kriebel  paper, 
although  they  are  also  found  in  the  Kranzer  and  Keller  report.  However, 

Fig.  1 differs  from  the  corresponding  graph  in  Kranzer  and  Keller.  In  the 
Kranzer  and  Keller  paper  the  maximum  value  of  A is  calculated  incorrectly, 
and  the  labels  on  the  graphs  do  not  indicate  correctly  what  is  actually 
being  plotted. 

The  function  n(r,t)  given  by  Eq.  (4)  could  be  dissected  into  several 
portions.  First  is  the  term  IQR^/pg^,  which  is  a constant.  It  is 
obvious  that  the  wave  height  should  be  proportional  to  the  magnitude  of  the 
initial  impulse.  It  is  also  proportional  to  the  square  root  of  the  effective 
radius  of  the  initial  impulse,  but  it  should  be  remembered  that  in  the  far 
field  the  initial  impulse  will  be  seen  as  a point  source.  Therefore,  as  r 
increases,  the  initial  Impulse  is  seen  to  be  more  and  more  like  a point 
source.  n(r,t)  is  seen  to  be  inversely  proportional  to  p and  g. 

Certainly  increasing  either  p or  g would  cause  an  increase  in  the  force 
against  the  upward  displacement  of  the  medium. 
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Secondly,  consider  the  inverse  proportionality  to  the  radius  r.  The 
amount  of  energy  in  the  entire  system  is  fixed  by  the  energy  contained  in 
the  initial  impulse.  As  the  wave  moves  out  radially,  this  energy  is  spread 
. out  over  greater  areas.  Its  intensity  diminishes,  and  so  the  wave  height 

decreases  accordingly. 

The  varying  amplitude  factor  A is  defined  in  Eq.  (7).  It  depends 

on  the  ratio  r/(gh)^2t,  as  seen  by  comparing  Eq.  (8)  and  Eq.  (7).  If 

r > (gh)ly,2t,  then  A is  seen  to  be  zero.  This  follows  from  the  fact 
that  the  maximum  group  velocity  of  any  wave  is  (gh)^2.  A plot  of  |A| 
is  given  in  Fig.  1 . 

Inspection  of  Eq.  (7)  indicates  that  A is  dependent  on  the  initial 
distribution  of  impulse  T(o/h),  the  effective  radius  R,  and  the  depth 
h.  All  three  must  be  specified  for  a particular  case.  Figure  1 plots  |A| 
versus 

1 . M1/?  t 

<P  r 

for  a parabolic  impulse 

= H ” 1 ( £ )21  for  r < 21/2  R 

=0  for  r > 21/2  R 

So  at  r = 0, 

ijrl  . , 

0 

and  at 

r = (2)1/2  R,  4^-  = 0 
lo 

The  effective  radius  R is  taken  as  (2/7 ) h and  h is  taken  as  5 
kilometers.  These  were  the  values  used  by  Kranzer  and  Keller  for  their 
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examples.  Within  their  paper  there  are  three  graphs  of  n(r,t)  versus  r 
or  versus  t (Figs.  4 and  5)  for  the  above  specified  values  of  R,  h, 
and  initial  impulse,  where  !Q  = 5 * 10  dyne/sec/cm  . These  same  values 
were  used  initially  in  the  project  in  order  to  verify  that  the  computer 
program  was  working  properly.  By  comparing  the  graphical  display  of  the 
computer  output  to  the  Kranzer  and  Keller  plots  it  was  found  that  although 
the  graphs  had  very  similar  amplitudes  and  wavelengths,  there  were  some 
discrepancies.  This  may  be  due  to  the  missing  a in  their  Eq.  (2.5),  as 
already  mentioned.  However,  their  Eq.  (2.11)  does  include  this  a.  The 
second  paragraph  of  Page  401  of  their  report  shows  that  they  wanted  to 
plot  A[r/(gh)  t]  on  the  abscissa.  The  maximum  value  of  |A|  given 
by  their  Fig.  3 and  Table  1 of  Amgx  = .68  is  incorrect  in  any  event. 

The  computer  program  gave  a maximum  value  of  A as  approximately  .96, 
and  this  was  the  same  value  reported  in  Kriebel's  paper.  Therefore,  the 
graphs  given  later  in  this  report,  which  have  the  same  initial  impulse, 
distances  or  times  that  Kranzer  and  Keller  indicated  in  their  Figs.  4 and  5 
are  to  be  taken  as  correct. 

The  remaining  term  in  Eq.  (4),  sin  2ir[(t/T)  - (r/X)],  is  the  portion 

of  the  wave  height  function  which  causes  the  more  rapid  oscillation  of  the 

wave  height.  The  varying  amplitude  A acts  as  an  envelope  for  this  more 

rapid  oscillation.  This  is  easily  seen  in  the  plots  of  n(r,t).  Both  the 

period  T and  the  wavelength  X are  dependent  on  the  dimensionless  ratio 
1/2 

r/(gh)  ' t.  T and  X are  not  fixed  constants,  so  the  argument  of  the 
sine  varies  in  a complex  manner.  When  r/[(gh)^  t]  = 1 , then  <j>(o)  = 1 
and 
★ 

This  refers  to  Fig.  3,  which  is  also  labeled  correctly  on  the  ordinate, 
but  the  caption  in  the  graph  indicates  its  inverse. 
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lim  1 t tanh  a ,1/2  . 1 (o)1/2  _ , 

o 1 0 7 ' 0 1 (TTnh  a)'/2  " 

So  when  o = 0,  then 

. - 2irh 

A = -V  00 

a 

and 

T ' (Sf  JL)1'2  * ‘ ' 

1 /2 

r [ ( gh ) t]  = 1 corresponds  to  the  outermost  part  of  the  wave  motion. 
Therefore,  the  wavelength  and  period  will  be  infinite  at  the  extremeties 
of  the  wave  pattern.  Both  T and  X will  decrease  as  r decreases 
for  a fixed  t,  or  when  t increases  for  a fixed  r.  It  is  to  be  expected 
that  T and  A should  increase  near  the  outer  portions  of  the  wave  motion. 
This  is  because  the  wave  can  be  considered  as  many  different  components 
which  emanate  from  the  origin  at  constant  but  at  many  different  speeds. 

Since  the  faster  portions  will  form  the  outermost  section  of  the  wave 
motion,  T and  X should  be  greatest  there.  This  follows  from  the  fact 
that  the  speed  of  a wave  varies  proportionally  with  its  wavelength  and 
period. 

From  the  above  discussion  it  is  clear  that  the  wave  pattern  expands 
as  it  propagates.  The  amount  of  this  spreading  out  is  proportional  to  t 
or  to  r.  For  example,  for  a given  t,  and  considering  the  wave  along  one 
ray,  the  distance  that  the  wave  pattern  occupies  at  large  r is  much  greater 
than  at  small  r. 

It  should  be  noted  that  A = 0 for  [r/(gh)1^t]  >1.  As  mentioned 
earlier,  r[(gh)^2  t]  = 1 determines  the  outermost  portion  of  the  wave 
motion,  as  at  distances  or  times  which  satisfy  [r/ (gh) 1 t]  > 1 , the 
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amplitude  of  the  wave  must  be  zero.  [r/(gh)^2  t]  > 1 corresponds  to 
p > 1. 

Looking  at  Fig.  1,  it  is  seen  that  |A|  goes  through  several  maxima. 
Had  the  scale  been  continued  to  1/<J>  °°,  there  would  have  been  an  infinite 

number  of  maxima.  The  distance  between  the  maxima  is  increasing,  and  this 
is  due  to  the  aforementioned  spreading  out  of  the  wave  pattern.  It  is  a 
linear  dependence  on  either  t or  r which  determines  this  distance 
between  the  maxima  because  <p  varies  linearly  with  r and  t,  as  seen 
by  Eq.  (8).  (It  should  be  noted  that  the  a which  satisfies  Eq.  (8)  is 
the  unique  non-negative  root  of  that  equation).  The  actual  value  for  the 
maximum  of  A will  be  discussed  later. 

Consider  now  the  initial  impulse  distribution.  In  this  project  the 
initial  impulse  distribution  is  parabolic  in  nature  and  is  mathematically 
described  by 


= 0 


for  r < 21/2  R 
for  r > 21/2  R 


Since  in  the  examples  used  h * 5 kilometers  and  R = ( 2/7) h , then 
R = 1.429  kilometers.  A plot  of  I(r)/IQ  is  shown  in  Fig.  4.  It  should 
be  noted  that  at  the  value  r = R,  the  impulse  has  the  value  [I(r)/I0]  = .50. 
The  effective  radius  is  defined  by 


\ |I  | RZ  - f r 1 1 (r ) | dr 
o 

Then  it|I0|R  Is  the  value  of  the  total  initial  Impulse,  and  IQ  is  the 
maximum  value  of  the  Impulse  distribution. 

In  order  to  compute  A from  Eq.  (7),  It  Is  necessary  to  know  the  value 
of  T(o/h)  where 
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1 ( F ) = / I(r)  J0(  IT  ) r dr 
o 

i.e.,  T(o/h)  is  the  zero  order  Hankel  transform  of  I(r). 
to  reduce  the  integral  as  follows: 

TlfM,  ;*I|rIJo(2r)rdr 

= ^ ‘ \ ( )2]  for  r < fZ  R 


It  is  possible 


for  r > R 


so  integral  (i)  becomes 


= / [1  - I ( £ )23  J0  ( T ) r dr  r 0 dr 


.2  ^ R 


'7/  '.<*&><*&>*<*) 


1 R , 

- -T  / r3  ( $¥■  ) dr 
2R2  0 0 h 


Recall  the  relationships 
/ x J (x)  dx  = x J,(x) 


/ xm  JQ(x)  dx  = xm  J^x)  + (m-1)  xm_1  JQ(x) 


(m-1 ) 2 / xm  2 JQ(x)  dx 


Equation  (ii)  becomes 


°o(x)  3 I Jl(x)  • J2(x) 


Then, 


I 
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Independent  of  the  shape  of  the  initial  Impulse  distribution  is  the 
actual  amount  of  energy  that  is  Imparted  to  the  medium.  The  amplitude  A 
will  asymptotically  approach  zero,  as  r/(gh)^  t approaches  zero,  if  the 
amount  of  energy  imparted  is  finite.  The  total  energy  Q is  given  by  the 
expression 

Q = l f°  C F T ( K )32  tanh  0 h d(oh)  (10) 

^ o 1 1 

Of  course.  In  all  real  situations  Q will  be  finite.  Only  if  I ( r ) is 
discontinuous  will  Q be  Infinite. 

Continuing  with  the  discussion  of  the  parabolic  Impulse,  it  is  possible 
to  find  the  maximum  values  of  all  the  variables.  Figure  3 shows  that 

-*•  /l  Jo  for  o > 4 (A) 

and  for  o > 4,  Fig.  2 shows  that  <J>(o)  -*•  — - — 

2 Jo 
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\ 

! 


M 


Then 

<fr 

Since 

<*> 


then 

a 


a 

4 


or 


| 


(B) 


Substituting  relationships  (A)  and  (B)  into  Eq.  (7),  then 


A % a S2o  T(g  h)  _ (o)3/2  (R)3/2  T(oh) 

h I M h3^2  I R2 

o o 

but, 

. Hi  ill  t3  . . ati  ,3/2 

(h)3'2  ' 43'2  r3  h3'2  ' 4r2 


so  the  variable  h plays  no  part  In  determining  the  amplitude  A for 
a >_  4.  Therefore,  the  depth  of  the  water  is  important  only  in  the  leading 
portion  of  the  wave  envelope,  but  the  frequency  of  the  actual  wave  is 
always  dependent  on  h. 

Introducing  the  variable  s * o/h  then  for  a > 4, 

A = /l  (Rs)3/2  -Jkl 
‘o  R 

and  s has  already  been  shown  to  be  independent  of  h for  a > 4.  For  a 
parabolic  impulse,  it  has  already  been  proven  that 


(2)  J2(  r?  Rs) 
(Rs)* 
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I 


Then, 


A 


ft  (Rs)3/2  (2)  J2(  ft  Rs) 

(rT? 


2 ^ (Rs)'1/2  J2(  ft  Rs) 


7 2 

Kriebel  indicates  that  at  Rs ( = gRt  /4r  ) = 2.10/ ft,  A is  at  a maximum. 
At  Rs  = 1.91,  J2(^  Rs)  = .473  and  (Rs)"1/2  = /OTT/1.64.  Then 
Amax  = ( 2 • 828) ( 1 . 41 4) 1 /2  (.286)  = .962.  This  is  the  first  maximum  indi- 
cated in  Fig.  1 . 

At  Rs  = 1.91,  Ra/h  = 1.91,  and  at  the  transition  point  o = 4, 

R/h  = 1.91/4  = .48,  and  ft  R = .68h.  So  for  Rs  £1.91,  the  depth  h 
has  no  effect  on  the  amplitude  A.  Therefore,  if  ft  R £ .68h,  then  the 
peak  value  of  A is  independent  of  the  depth  of  the  medium.  ( ft  R is 
the  actual  radius  of  the  initial  impulse.) 

In  summary,  the  following  maximum  values  exist. 


A = 
max 

°max 

*max 

Tmax 

t 

max 

^ax 

Kranzer 


.962 

1.91  h/R  = 

OTR  = 2-3 

3.7  (R/g)1/2 

r) 

•962  lo  rr 

p r \l  g 

and  Keller  also  mention  other  distributions  than  parabolic  ones. 


Their  Table  I has  been  reproduced  below.  Listed  are  the  maximum  values  of 
some  of  the  variables  for  five  different  distributions  of  the  initial 
Impulse.  Table  II  and  Fig.  5 contain  some  values  for  the  same  distributions 
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TABLE  I.  CHARACTERISTICS  OF  FIVE  DIFFERENT  DISTRIBUTIONS  OF 
INITIAL  IMPULSE  (From  Kranzer  & Keller) 


Hr) 

I, 

1,  r<K 
0,  r>K 

NG)’}  's!'" 

0,  r>2*« 

hG)T ' 

t-  (./*/> 

<'l  mu 

1.1 

0.96 

0.49 

0.42 

0.70 

<T« 

0.33 

0.30 

0.35 

0.42 

0.38 

K 

2.8 

23 

3.1 

4.4 

3.6 

r..„ 

(*/*)» 

4.1 

3.7 

4.3 

5.3 

4.8 

Sccomlary 

maxima 

Infinitely  many;  practically 
constant  amplitude 

Infinitely  many;  ampli- 
tudes decrease  like  r* 

.Vonc 

None 

None 

Q/rRfT'If 

•0 

0.48 

0.18 

0.14 

0.31 

IU) 

(»—/<-') 

/. 

> R 

r.fp.\ 

2 

rsjiPwi 

j ftly-lt./ll  . 

s 

/* 

2[1  + J/Pj»3I 

TABLE  II  (From  Kriebel) 


I/I 

O 

A 

m 

Eq.  (B.45) 

m 

A / Vq 
m 

1(0/ 1 

o 

Eq.  (B.l) 

I o 
o 

Table  R-l 

i 

1.16 

OC 

0 

- J, (Ks) 
s 

to 

2 

0.96 

0.48 

1.39 

~J  J2/2Rs 

/2 

3 

0.49 

0.18 

1.15 

r2  - £5. 
2~e  /2 

2/2 

4 

0.42 

0.14 

1.12 

_a 

R2  / R2s2\  2 

r i1  * — ) 

00 

5 

3.76 

0.31 

1 .37 

y.  ■ (?)' 

/ir 
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as  given  by  Kriebel.  Included  in  Fig.  5 is  a graphical  display  of  the  dis- 
tributions. The  parabolic  impulse  has  been  labeled  as  number  2.  (Note 
that  Vmax  refers  to  the  group  velocity  in  Table  II.) 

Among  the  distributions  listed  is  a discontinuous  distribution  given 
by 

« 1 r < R 

xo 

= 0 r > R 

As  a result  of  the  discontinuity,  Q is  infinite.  This  causes  infinitely 

many  maxima  to  be  generated  which  do  not  approach  a zero  magnitude,  and  the 

distance  between  maxima  decreases  as  t increases.  The  other  four  distribu- 
tions are  all  continuous  and  so  contain  a finite  amount  of  energy.  Kranzer 
and  Keller  point  out  that  "The  last  three  have  no  secondary  maxima  while 
the  first  has  infinitely  many,  the  magnitudes  of  which  decrease  like  t 
raised  to  the  minus  2 power  (for  fixed  r).  The  existence  of  secondary 
maxima  and  their  rate  of  decrease  depend  upon  the  smoothness  of  I(r).  If 
I ( r ) Is  discontinuous  infinitely  many  maxima  of  practically  constant  ampli- 
tude will  occur.  If  I ( r ) is  continuous  but  its  derivative  discontinuous 
for  some  r f 0,  infinitely  many  maxima  will  occur  with  amplitudes  decreas- 
ing like  t raised  to  the  minus  2 power  (for  fixed  r).  If  I(r)  and  all 

its  derivatives  are  continuous  there  are  usually  only  a finite  number  of 
maxima,  and  If  I (r ) is  a decreasing  function  only  one  maximum  can  be 
expected."  (from  p.  401). 

In  the  Kriebel  paper,  the  case  of  very  deep  water  (h  <*>)  is  also 
considered.  By  examining  this  situation  the  effect  of  changing  the  effective 
radius  R and  the  total  Impulse  and  total  energy  can  be  found. 

Previously,  in  Eq.  (10),  the  quantity  Q was  defined  as  the  total 


energy  imparted  to  the  medium.  When  considering  a parabolic  distribution, 
with  h -*•  «,  and  using  the  relationship  of  Eq.  (9),  then. 
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Q - J22(/?  Rs)  d(Rs) 

-2JL.  = 4 / , 

" RI  2 o (Rs)2 


. 16 

B -tt—  * .480 


Then, 


.480  tt  R I 


The  total  impulse  was  also  previously  defined  and  is  given  by 


J - * R‘  i0 

Equation  (4)  defines  n(r,t)  as 
I R1/2 

n(r,t)  % °-  m — A sin  2ir  ( j-  - £ ) 


P 9 r 


It  has  been  shown  that  Amflx  = .962,  and  the  maximum  value  of 


sin  2ir  [(t/T)  - (r/X)]  = 1.  Then, 


.962  IQ  R1 


Therefore,  in  terms  of  J and  Q, 

n - .962  j xnr 

%,ax  V^Vr  V 9 

_ .785  ( ~Q 
^max  " r | p g 


.48  p g 


-3/2 

Then  for  a given  total  impulse  J,  nmax  is  proportional  to  R . For 
a given  Q,  nmax  Is  Inversely  proportional  to  r and  is  independent  of 
R.  (Kriebel  points  out  that  the  time  of  the  arrival  of  this  maximum  wave 
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height  is  t = /(7.64/gR)  r,  which  is  dependent  on  R.) 

From  Fig.  5 it  is  seen  that  as  I ( r ) becomes  more  peaked,  then  the 
value  of  Amax,  and  therefore  nmax.  decrease.  For  example,  I/IQ  of 
curve  4 is  more  peaked  than  curve  2,  and  Amav  is  almost  half  as  great 
for  curve  4.  It  must  be  remembered  that  all  of  these  curves  and  all  com- 
parisons are  made  with  the  same  IQ  and  R.  Ans  since  both  J and  Q 
are  determined  by  only  these  two  variables,  then  the  total  impulse  and 
total  energy  in  all  of  these  cases  is  the  same. 

DISCUSSION  OF  THE  COMPUTER  PROGRAM 

Using  a CDC  7600  computer  a program  caled  WWAVES  was  written  to 

evaluate  the  function  n(r,t).  Basically,  Eq.  (4)  was  solved  using  the 

values  calculated  in  Eqs.  (5)  - (8). 

Kranzer  and  Keller  suggested  that  a value  of  $ be  selected,  and  then 
the  corresponding  values  of  either  r or  t calculated  from  Eq.  (8). 

This  was  the  first  method  in  approaching  this  project;  however,  this  led  to 
equal  increments  in  $ and  varying  increments  in  r or  t.  Since  it  was 

desired  to  later  also  superimpose  the  waves  resulting  from  several  initial 
impulses,  it  was  necessary  to  have  equal  increments  in  r or  t. 

To  find  equal  increments  in  r (or  t),  Eq.  (8)  was  solved  for  o for 

a given  t (or  r).  Equation  (8)  is  a transcendental  equation  in  a and 
not  readily  solved.  Figure  2 graphs  Eq.  (8).  Consider  now  a fixed  t and 

r varying,  r will  be  set  at  some  lower  limit,  which  satisfies  r « R. 

In  order  now  to  solve  for  o,  Newton's  method  was  first  applied.  But  due 
to  the  flatness  of  Eq.  (8)  as  seen  in  Fig.  2,  this  method  of  root  finding 
was  very  slowly  converging  in  some  regions  and  even  diverging  in  some  areas. 
Next,  the  Method  of  False  Position  was  used  with  greater  success.  A schematic 


diagram  of  this  portion  of  the  program  is  found  in  Fig.  6.  It  is  based  on 
the  equation 


19 


X(n+1 ) 


‘(n) 


, lf(x(n)> 1 

lf(x(n-l)) I + lf(x(n)>l 


'(n) 


) 


This  schematic  is  taken  from  page  10  of  An  Introduction  to  Determination 
of  Roots  of  Equations,  by  C.  D.  Mote,  Jr.  Quoting  from  page  8 of  these 
notes,  "The  method  (of  False  Positon)  is  used  to  determine  a real  root 
xi  of  the  equation  f(x)  = 0 when  values  of  the  function  f(x)  are 
known  at  x„  = x-  - and  x,„  ,\  = x4  + . The  root  is  estimated  by 

linearly  interpolating  between  the  values  of  the  function  (with  opposite 
signs)  at  successive  trial  positions.:  Graphically: 


The  disadvantage  of  using  the  method  of  false  position  is  that  the 
root  must  be  between  two  known  values.  Therefore,  values  of  a were  pre- 
selected between  which  the  root  value  of  a was  definitely  known  to  lie. 
This  restriction  and  the  flatness  of  the  equation  caused  slow  convergence. 
The  method  of  false  position  is  an  iterative  one,  and  the  test  for  con- 
vergence was  set  to  a pre-selected  tolerance.  For  all  of  the  output  in 
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this  report  the  error  criteria  is  .01.  Greater  accuracy  is  easily  available, 
but  involves  more  computer  time. 

Once  a was  found,  it  was  substituted  into  Eqs.  (5),  (6),  and  (7), 
and  these  resulting  values  allowed  n(r,t)  to  be  calculated.  Then  the 
value  of  the  varying  variable  was  incremented  and  the  cycle  started  over 
again. 

Next  the  maximum  and  minimum  values  of  n(r,t)  and  either  r or  t 
were  found.  The  general  form  of  the  printed  program  can  be  found  later  in 
this  report. 

Now  it  is  possible  to  use  the  Graphical  Display  System  (GDS)  which  is 
in  the  computer  library  to  plot  n(r,t).  The  necessary  parameters  for  the 
plot  are  set  and  the  graph  is  prir*ed  on  the  line  printer.  The  GDS  file  is 
then  disposed  to  microfiche,  whic  jives  a fine  line  plot.  Only  microfiche 
copies  of  the  output  are  found  within  this  report.  The  subroutine  available 
in  the  GDS  library  is  entitled  PFLILI.  It  utilizes  a parabolic  fairing 
technique  to  draw  a curve  through  a set  of  points.  Therefore,  the  plot  of 
n(r,t)  is  a very  smooth  curve  and  not  a series  of  straight  lines.  Also 
included  in  the  plot  is  the  wave  envelope,  indicated  by  dotted  lines. 

This  section  of  the  report  includes  several  cases  of  initial  impulse 
distribution.  The  first  case  includes  the  computer  output  partially  show- 
ing the  digital  results.  The  corresponding  plot  graphically  displays  the 
wave  neight  versus  time  for  three  initial  impulses.  The  resulting  waves 
from  the  initial  impulses  have  been  superimposed.  It  should  be  noted  that 
since  the  size  of  the  graph  on  the  vertical  axis  is  constant  for  all  the 
graphs,  it  is  possible  to  compare  wave  heights  between  different  cases  only 
by  comparing  their  actual  value  and  not  their  height  on  the  graph.  Since 
it  is  possible  to  superimpose  several  waves,  a study  of  the  waves  resulting 
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from  actual  impulses,  rather  than  idealized  parabolic  ones,  can  be  under- 
taken . 

The  next  six  plots  display  wave  height  versus  time  for  one  particular 
initial  impulse  while  the  distance  from  the  impulse  has  been  varied.  Note 
that  the  wave  envelope  spreads  out  as  the  distance  from  the  impulse  changes 
from  20  to  30  kilometers. 

The  next  set  of  graphs  are  of  wave  height  versus  distance  with  time 
varying  from  350  to  500  seconds.  As  the  time  from  the  initial  impulse  is 
increased,  the  wave  envelope  can  be  seen  to  be  spreading. 

Finally,  the  last  two  pages  contain  a comparison  between  empirical 
data,  and  information  generated  by  the  computer  program.  Full  information 
for  the  actual  tests  was  not  available,  so  it  is  impossible  to  fully  compare 
the  theoretical  with  the  actual  results.  In  the  second  of  the  two  experi- 
ments cited,  the  depth  of  the  water  was  not  given,  but  was  assumed  to  be 
130  feet  (an  approximation  for  deep  water).  This  approximation  may  have 
caused  the  discrepancy  in  the  frequencies.  It  should  be  noted  that  in  the 
first  graph  presented  the  wave  does  not  touch  the  wave  envelope.  This  is 
due  to  the  low  number  of  data  points  that  the  computer  was  able  to  generate. 
It  is  felt  that  what  is  really  important  to  gain  from  this  graph  is  the 
value  of  the  maximum  of  the  wave  envelope.  This  also  applies  to  the  second 
graph,  although  more  resolution  could  have  been  obtained  had  the  program 
been  instructed  to  generate  more  data  points.  It  was  felt  that  the  time  and 
expense  involved  was  not  necessary.  In  general,  it  can  be  concluded  that 
there  was  fairly  good  correlation  between  the  empirical  data  and  the  wave 
shapes  predicted  by  Kranzer  and  Keller's  theory. 
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versus  1/<J>  FOR  A PARABOLIC  IMPULSE,  R = (2/7)h,  h = 5 KM,  and  I = 5 x 10'  dyne  sec/cm1 
(from  Kriebel) 
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FIG.  6.  FALSE  POSITION  COMPUTER  FLOW  DIAGRAM 
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FIG.  17. 
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The  figure  below  is  the  wave  envelope  and  waves  generated  by  an  explosion 
of  TNT  on  the  ocean's  surface.  It  was  calculated  using  Kranzer  and 
Keller's  Theory,  but  is  to  be  compared  to  an  actual  test  run  on  Mono 
Lake  in  1965.  The  radius  of  the  charge  was  34  inches  (weighing  approxi- 
mately 5 tons),  and  the  depth  of  the  water  was  130  feet  (.04537  kilometers) 
Equation  B.47  of  Kriebel  was  used  to  calculate  the  effective  radius  of  the 
explosion  and  the  initial  impulse  impared  to  the  water.  Equation  B.47  is: 

R = (2)1/3  x .17  x a = (2)1/3  x 11  x RQ 

where  a = 64.5  x RQ  for  TNT;  RQ  is  the  charge  radius;  R is  the  effec- 
tive radius 

(2),/3  . pa 

I = - 

0 

where  p.  is  atmospheric  pressure;  c,  is  the  speed  of  sound  in  air; 

a a 

a = 64.5  x Rq  for  TNT,  and  IQ  imparted  to  the  water  is  ten  times  IQ 
imparted  to  the  waves. 

In  the  actual  test,  nmax(r)  = 1800  ft2,  so  at  200  feet, 

\ax  = 2’^  meters-  Kranzer  and  Keller's  result  is  fairly  close  to  this 
result,  with  most  of  the  error  probably  coming  into  the  results  through 
the  approximation  for  the  initial  impulse. 
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FIG.  18. 
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Comparison  of  theoretical  and  experimental  data  is  shown  below.  The  upper 
graph  is  the  wave  envelope  measured  200  ft.  away  from  a 125  pound  charge 
in  1963  (Pinkston  1966).  It  is  taken  from  Kriebel , Fig.  A-2.  The  lower 
graph  is  generated  using  Kranzer  and  Keller's  theory,  with  the  initial 
impulse  calculated  using  Eq.  B.47  of  Kriebel.  (Note  that  Eq.  B.47  gives 
the  impulse  imparted  to  the  waves,  but  the  impulse  imparted  to  the  water 
is  ten  times  this  amount.) 

R = 0.675  ft 
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